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Abstract 

Density fluctuations of expanding nuclear matter are studied within a mean-field model in which 

fluctuations are generated by an external stochastic field. The time evolution of the system is 

studied in a kinetic-theory approach. In this model fluctuations develop about a mean one-body 

'■^ , phase-space density corresponding to a hydrodynamic motion that describes a slow expansion of the 

\ system. A fluctuation-dissipation relation suitable for a uniformly expanding medium is obtained 

' and used to constrain the strength of the stochastic field. The coupling between the kinetics 

of fluctuations and the hydrodynamic expansion is analyzed, and the distribution of the liquid 

domains in the spinodal decomposition of this expanding nuclear matter is derived. It is found 

, that the formation of the domains can be envisaged as a stationary process. Comparison of the 

related distribution of the fragment size with experimental data on the nuclear multifragmentation 

CnI ' is quite satisfactory. 
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I. INTRODUCTION 



In a previous paper [lj we have studied a model in which the excited nuclear system that 
is formed after the collision of two heavy ions at intermediate energy is described as hot 
unstable nuclear matter. In this system the density fluctuations can lead to a liquid-gas 
phase transition ( spinodal decomposition ). The size of the liquid domains compares well 
with the experimentally observed yield of fragments. In the model of Ref. |l| the density is 
assumed to fluctuate about a stationary mean value qq. In order to take into account the 
expansion of the complex system formed in the reaction, here we generalize that approach 
and consider the more realistic situation in which the density go can be time-dependent: 
Qo — Qo{t)- We assume that the mean density Qo(t) changes slowly with time, compared to 
the times characterizing the behavior of fluctuations, and obtain a hydrodynamic description 
of the time dependence of Qo(t). In order to describe fluctuations about this time-dependent 
average behavior, we assume, like in Ref. jl| , that a stochastic self-consistent mean field in- 
duces fluctuations of the density about its average value. The evolution of these fluctuations 
is accordingly described by means of a kinetic equation of mean-field in which a stochastic 
( Langevin ) term is included. 

Semiclassical kinetic equations for the one-body phase-space density have often been used 
for studying the complex processes occurring in heavy ion collisions 0,0,0). However, in 
their original formulation, these equations account for the time evolution of the average one- 
body phase-space density only, and cannot describe phenomena such as multifragmentation 
in which fluctuations of the phase-space density about its mean value are believed toplay 
an essential role ( for a review on nuclear multifragmentation see, e.g., Refs. 0, E| 0, 
). The Boltzmann-Langevin equation of Refs. [2, El 11| instead, incorporates also 



a stochastic term into the kinetic equation and includes fluctuations in the evolution of 



3QUatlO] 

the phase-space density. In Refs. (§, EH LL1( the diffusion coefficient of the Langevin 



fluctuating ) term is ultimately related to the amplitude of the nucleon-nucleon scattering. 
That relation is a particular case of the fluctuation-dissipation theorem. More recently, a 
different method to take into account fluctuations has been proposed in Ref. In that 



paper the statistical fluctuations of the one-body phase-space density are directly introduced 
by assuming local thermodynamic equilibrium. A basic assumption, shared by all works on 
this subject, concerns the white-noise nature of the stochastic term. In Ref. |l| we have 
shown that this assumption is compatible with the the fluctuation-dissipation theorem only 
in particular situations since this theorem provides a kind of self-consistency relation for the 
stochastic field. That result is generalized here to the case of a slowly expanding system. 

Thus, we extend the approach of Ref. [l| by considering fluctuations about the average 
time-dependent density Qo{t) of nuclear matter during a hydrodynamic expansion. The 
evolution of the fluctuations is described by a kinetic equation that in the absence of the 
stochastic field has a solution corresponding to the one-body phase-space density of the 
hydrodynamic motion. The stochastic field is treated as a linear perturbation about the 
zero-order hydrodynamic solution. The diffusion coefficient of the stochastic field is self- 
consistently determined by means of a fluctuation-dissipation relation, which is proved to 
be valid for a uniformly expanding system. We are mainly interested in unstable situations, 
where the fluctuations increase with time until they cause the decomposition of the system. 
The growing of fluctuations is essentially dominated by the unstable mean field. Thus 
we focus our attention on the behavior of the mean field and neglect the collision term 
in the kinetic equation. Collisions would mainly add a damping to the growth rate of 
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the fluctuations and could not change the main results of our calculations, at least at a 
qualitative level. 

In Sec. |n] we outline a procedure to determine the structure of liquid domains formed 
within the system during its spinodal decomposition. This formalism is a generalization 
of that introduced in Ref. [lj and it gives a new insight into the interplay between the 
hydrodynamic motion and the kinetics of fluctuations. The present approach allows us to 
determine the time behavior of the distribution of the liquid domains, that depends on the 
collective expansion energy. Moreover, our results for the expanding system reproduce the 
main results of Ref. 

According to Ref. [13[ the fragmentation phenomenon observed in heavy ion collisions 
could be ascribed to a spinodal decomposition of the bulk of nuclei. Following that sug- 
gestion, we identify the pattern of domains formed in the decomposition of our expandig 
nuclear matter with the fragments experimentally observed in multifragmentation reactions 
and relate the distribution of the liquid domains to that of the fragments. Thus, in Sec. 11111 
we compare the results of our calculations with suitable experimental data generated by 
processes for which there are indications that fragmentation could ensue from the expansion 
and cooling of the emitting source, following either an initial compression 0, ^| or an 
initial heating rJ n|. Finally in Sec. IV a brief summary and conclusions are given, while 



in the Appendix the statistical properties of expanding nuclear matter are briefly discussed. 



II. FORMALISM 

A. Kinetic equation 

We want to study the spinodal instabilities of nuclear matter when it is brought into the 
spinodal zone of the phase diagram. For this purpose we model the system formed after the 
collision by a sphere of homogeneous nuclear matter that is slowly expanding and assume 
that the radius of this sphere is much larger than the characteristic bulk lengths. Then we 
can study bulk properties of nuclear matter and confine our attention to the core of the 
sphere by neglecting both surface and finite-size effects. We further suppose that the core 
of the excited system is made of homogeneous nuclear matter with density and temperature 
that are functions of time only. Hence the continuity equation 

^ + eo (f)V.u(r,t) = 0, (1) 

requires that the divergence of the hydrodynamic velocity field V • u(r, t) also depends only 
on time. Choosing the form 

uM) = -rg| (2) 



for the hydrodinamic velocity, Eq. gives 

Qo(t) _ l 3 (t) 
00 (0) 7 3 (0) 

for the time evolution of the density. 



(3) 
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The function j(t) is determined by the second Euler equation 

JU(r,t) + (u(r,t)- V)u(r,t) = 0, (4) 

which is fulfilled with 

^ = 1 (5) 
7(0) (1 + at) ' 1 J 

where a = |t(0) /t(0) | is the initial expansion rate. 

Since the velocity field is simply proportional to r, both shear and bulk viscosities do 
not give rise to any effect. Then, we can assume that the nuclear system moves along an 
isoentrope during the expansion. From a microscopic point of view, we shall treat nuclear 
matter within a mean-field approximation. The one-body phase-space density appropriate 
to the present hydrodynamic motion is 

io(p, r, t) - e/9(t) ( (p _ mu(r)t))2/(2m) _- w ) + 1 ' ( 6 ) 

where j3(t) = 1/T(t) is the inverse temperature ( we use units such that h = c = k B = 1 ), 
and the effective chemical potential ft(t) is measured with respect to the uniform mean field 
U(@o{t))- Within this scheme the relations among temperature, effective chemical potential 
and density along the isoentrope are the same as for a Fermi gas nana 

T(t) _ jj(t) _ / g0 (tK2/3 

r(o) mo) U(oy ' 1 ) 

Thus, in our model, the average trajectory of the system in phase space is determined by 
the hydrodynamic expansion. Following Ref. |]J, we introduce now a stochastic scalar field 
that generates fluctuations about the average trajectory and calculate the density-density 
response function of the system in a self-consistent mean-field approximation. In order to 
derive compact analytic expressions, we use the linearized Vlasov equation for calculating 
the response function. Thus we write the one-body distribution function when the stochastic 
field is present as 

/(p, r, t) = / (p, r, t) + 5f(p, r, t) 
and assume that at all relevant times 5f << Jq. Then 5f obeys the equation 

ot m 

+—fo(p, r, t) ( P - u) • V / A(r - r', t)5g(v', t)dv' 

= - 7F fo(p,v,t)(V-u)-Vv(r,t), (8) 

where f(r,t) is the external scalar field and Sg(r, t) is the density fluctuation 

5g(r,t) = g(r,t) - g (t) = J 5f(p,r,t)dp. 
The effective interaction ^4(r — r', t) is given by the functional derivative of the mean field 
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evaluated at the actual density of the expanding system. Here we use a schematic Skyrme- 
like, density-dependent, finite-range effective interaction that has been introduced in Ref. 
jiol ]. In momentum space it reads 

A k = (a— + (a + y-L-ffie-* k * /2 , (10) 

Qeq Q eq 

with 

a = -356.8 MeV, b = 303.9 MeV, a = - . 

6 

These values reproduce the binding energy ( 15.75 MeV ) of nuclear matter at saturation 
( Qeq — 0.16 fm -3 ) and give an incompressibility modulus of 201 MeV. 

The width of the Gaussian in Eq. (fTUj) has been chosen in order to reproduce the surface- 
energy term as prescribed in Ref. . For small k this interaction reduces to that employed 
in Ref. £0] and the parameter c in Eq. (jlOj) is obviously related to the parameter d of Ref. 

Q. 

In the physical situations considered here the values of temperature are sufficiently small, 
with respect to the Fermi energy, so that the Pauli principle is still operating. For T 
sufficiently low with respect to jl, the derivative df (p,r,t)/dfl in Eq. (JHJ) is appreciably 
different from zero only in a small domain of the nucleon velocity |p|/ra about the Fermi 
velocity Vp. Moreover, the expansion velocity |u| must be smaller than vf, otherwise our 
hydrodynamic approximation to the expansion would not be reasonable, so we can put 
(p/m — u) pa p/m in Eq. (jSJ). Neglecting the expansion velocity in Eq. (jSJ) simplifies 
calculations and, when taking the space Fourier transfom, leads to a closed equation for 
each Fourier coefficient. 

The integral equation satisfied by the response function in momentum space reads: 

D k (t,t')=Df\t-t',t l )+ f Df\t-t"X)A k {t")D k {t\t')dt". (11) 

J -to 

The origin of time is set at t = —to, with to much larger than the typical response times 
of the physical system. The non-interacting particle-hole propag ator D^\t-t',f) has the 
same expression as for a system at equilibrium, however in that case D% depends only on 
the time difference t — t' while here it gets a further dependence on t' through the time- 
dependence of the thermodynamic quantities that determine the local-equilibrium state. For 
symmetry reasons both propagators in (fTTj) depend only on the magnitude k of the wawe 
vector. Equation (fTTj) gives the response to a scalar field of a homogeneous system with 
thermodynamic properties that change in time according to an isoentropic expansion. For a 
vector field, instead, the hydrodynamic velocity field u(r, t) could not always be neglected, 
and spatial uniformity could be lacking. However, even for the density-density response, 
in our approach the kinetic evolution of fluctuations is still coupled to the hydrodynamic 
motion. 

Spinodal instabilities in expanding nuclear matter have been previously investigated in 
Ref. ji^], where the authors focused their attention mainly on the growth of the most 
unstable modes. In that paper, time-averaged values were taken for the quantities that 
change on the hydrodynamic time scale, so that the coupling between hydrodynamic motion 
and kinetic evolution of fluctuations was not fully taken into account. 
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B. Distribution of fluctuations 



According to the approach of Ref. we assume that the time scale of terms not included 
in the mean-field approximation is shorter than the characteristic times of the mean-field 
dynamics and treat all the more complicated processes like thermodynamic fluctuations, 
quantum effects and short-range correlations, as an extra stochastic field similar to the 
random force of the Langevin equation for Brownian motion. Moreover we assume that this 
additional field is a Gaussian white noise with vanishing mean so that the time-evolution 
of the density is a Markovian process. Like in Ref. pj, we shall determine the conditions 
for which the white-noise assumption can be considered valid also in this time-dependent 
environment. 

The additional stochastic mean field will induce density fluctuations with respect to the 
uniform density of the expanding system. To be more specific, we assume that at the time 
t = in the system is present a density fluctuation 5g(r, t = 0). We can imagine that 
this initial fluctuation is the result of a fictitious external force that has been acting from 
the time —t until time t = ( see for example Ref. [23], Sec. 15 I ) and assume that 
density fluctuations are absent before time —to- In linear approximation for the stochastic 
mean-field, the Fourier coefficients of 5g(r, t) for t > are given by 

5ek ® = /n^w* f ^ dt ' + f *WO<Wk(0 • (12) 

]_ to D k (0,t')dt' J -t Jo 

In the second integral, Bk(t')dWk(t') gives the contribution of the stochastic field in the 
interval dt'. The real and imaginary parts of the Fourier coefficients W^t') are indipendent 
components of a multivariate Wiener process [24]. Since the stochastic field is real B^(t) = 
B_ k (t) aadW2(t) = W_ k (t). 

In Ref. |l| we have shown that the white-noise hypothesis for the stochastic field can 
be retained for values of temperature and density sufficiently close to the borders of the 
spinodal region. Moreover, in the physical situations considered both in Ref. £0] and here, 
the strength of particle-hole excitations having energy larger than kvp can be neglected. 
The imaginary part of the Fourier-transformed response function Dk(uj,t') displays a sharp 
peak dominating the particle-hole background at a value of uj <C kvp. As a function of the 
complex tu-variable, the response function has a pole on the imaginary axis, at a distance 
from the origin that is smaller than the values of kvp. The position of this pole determines 
the time scale characteristic of the response function. This time is much longer than the 
times l/(kvp) characteristic of the non interacting propagator Df\t — t',t'). Thus, when 
taking the inverse Fourier transform, we can neglect the contributions from the cut along 
the real-cu axis and keep only that of the isolated pole. Moreover, here we are interested in 
density fluctuations that either relax towards their equilibrium values, or grow indefinitely 
in the unstable case. Thus we look for solutions of Eq. with values of (t — t') of the 
same order of magnitude as the damping or growth time scales. In this case we can neglect 
the inhomogeneous term in Eq. 

When our combined hydrodynamic and kinetic approach is justified, the time typical of 
the hydrodynamic expansion \^{t) / ^{t)\ is also longer than l/(kvp), as a consequence the 
integral equation (fTTj) is reduced to a simpler differential equation. Since the propagator 
Dfr(t — t",t") in the integral in (Jll|) depends on the second argument only through the 
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slowly changing hydrodynamic quantities, we expand it as 

D®(t - t"X) ~ D^\t - t",t) + (t- t")^Df\t - t",t) , (13) 
and neglect higher order terms. Thus we only need the two quantities 

DSV = o,t) = -»g9 

and 

OUO UK 

where the function 



e-Pi* + 1 

is constant during the isoentropic expansion, because of the relation (jJJ). 

With these replacements, the integral equation (|11|) turns into the following differential 
equation 



- —F(Pfi)A k (t)-D k (t,t') = (l + J^A k (t) - - —F(^)—^jD k (t,t') . (14) 
The solution of this equation can be put in the form 

D k {t^ = C k {t')^e^'^\ (15) 

where the function 

tt i an m f^pr^ + Ak ^ ~ Mt)) 

Lk[t) ~^FW)^ Mt) k {b) 

is the damping or growing rate ( depending on its sign ) of the density fluctuations. This 
function extends the result given in Eq. (2.16) of Ref. JlJ to a slowly changing system. The 
quantity f(t) in this equation is the free-energy density. The sign of T k (t) is negative for 
#7(0 , 
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\t > 0, corresponding to a system at equilibrium, while it is positive for 

d 2 f(t) 



del 



|r < -(A k (t) -A (t)) < 



when the system is inside the spinodal region. 

In order to derive Eq. (|16|) we have used the following relation 
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The function Ck{t') in Eq. ()15|) is still to be determined, for the moment we assume that 
it varies according to the hydrodynamic time-scale and shall verify a posteriori the validity 
of this conjecture. 

If the approximate solution (|T5j) for D k (t,t') is replaced into Eq. (|12|). the following 
expression for the density fluctuations 5g k (t) is obtained: 

op k {t)= 5g k {t = 0)———eJo k[ - > 

A.h{t) 

1 r* r* 



[ e ti v ^ dt ''A k (t')C k (t')B k (t') dW k (t') . (18) 
Jo 



A k (t) Jo 

This expression can be viewed as a solution of the stochastic differential equation 

dz k {t) = T k {t)z k {t)dt + A k {t)C k {t)B k {t)dW k {t) , (19) 

with z k {t) = 5g k (t)A k (t). The corresponding Fokker-Planck equation for the probability 
distribution P[z k (t)] reads 

+iA(() 2 C t (i) 2 |-B k (*)rg|^yPfe(*)]. (20) 

The stochastic part of the mean field is completely determined once the coefficients B k (t) 
are known. In order to obtain information about these coefficients we concentrate our atten- 
tion on the correlations of density fluctuations at local equilibrium. Due to the independence 
of the components of the multivariate Wiener process W k (t), only the correlation terms with 
k' = — k survive. The equation for the equilibrium fluctuations can be obtained from Eq. 
(112)) by shifting the initial time t = to —to, when the stochastic field has been turned on, 
without including any particular condition at later times. Then, the correlations are given 
by 

< Sg k {t)5g^{t r ) > ' 



A k (t)A k (t') 



J -t 

where the brackets denote ensemble averaging. 

The aim of our approach is that of determining the Fourier coefficients B k (t) in Eq. ()12j) 
as functions of £> (^) an d T(t) for the system at local equilibrium, and of extending afterwards 
the expressions so found to non-equilibrium cases. Such a procedure is usually followed when 
trea ting instabilities by making use of the fluctuation-dissipation theorem, see e.g. Refs. 



|25l 26]. As a part of this program, we have to verify the self-consistency of the withe- 
noise hypothesis for the stochastic field, i.e. the coefficients B k (t) must be local functions 
of time. This step can be accomplished by exploiting the fluctuation-dissipation relation 
obtained in the Appendix within the same approximations used in deriving Eq. ()11|) . The 
fluctuation-dissipation relation involves the Fourier transforms of the density correlator and 
of the response function with respect to the difference of times, whereas both the average 
value of the density correlations and the response function refer to a state that changes 



with time according to the hydro dynamic time scale \'j{t)/'y(t)\. The condition of local or 
instantaneous equilibrium for the hydrodynamic expansion requires that the kinetic time 
scales, l/|rfc(i)| and 1/vpk, are smaller than |7(t)/7(t)|. Thus, in Eq. (|2*Tj) the quantities 
Ak(t), Ck(t), \B^(t)\ and T k (t) are slow functions of time, compared to the exponentials, and 
can be approximated by the value that they take at the upper integration extreme. Within 
this approximation Eq. (|21j) gives ( choosing t > t' ) 

AJt) P r k (t')(t-t') 
< **(«) WO >= 4^.(0 W\ • (22) 

In Ref. |H we have shown that, if the classical limit uj/T <C 1 ( or \T}~(t)\/T <C 1 ) can 
be taken when evaluating terms on both sides of the fluctuation-dissipation relation, the 
assumption of a white-noise stochastic field can be retained. Since the relevant values of the 
wave vector k turn out to be such that the quantity kvF is of the same order of magnitude 
as T, the limit u/T <C 1 also implies u/kvp <C 1. In analogy with Ref. here we also 
assume that T > so that the result obtained there, concerning the validity of this 

inequality when temperature and density are in the proximity of the spinodal region, still 
applies. Then the fluctuation-dissipation relation of Eq. (jA.5j) can be put in the form 

^ < 5 6k (t)5g^(t') >= -T(t')D k (t, . (23) 

The variance of the density fluctuations of nuclear matter at local equilibrium can be 
determined by using Eq. (|A.4|) . Within the present approach ( see also Ref. Q ) it is given 
by 

Tit) 

al(t)U =< te ( t )^(0 >= nt) + _ Mt)) ■ (24) 

where we have introduced the abbreviation 



f" 



07, 

00 



JIT 



By exploiting this equation and Eq. (|23j) . we can obtain the following relation between the 
function Ck(t) and the diffusion coefficient |-Bk(^)| 

Cl(() = 2 wk- (25) 

and the additional relation 

\B k (t)\ 2 = -m 2 T(t)F(mMt)(A (t) - /"(*))- . (26) 

7T Kj 

These equations confirm that both |i?k(t)| and C k it) change in time according to the hydro- 
dynamic time-scale, as we have previoulsly assumed. 

Both the diffusion coefficients of Eq. (f2T)J) and the function of Eq. (|23|) have been derived 
for a system at local equilibrium using the subsidiary conditions j3 < l/II^I < I7/7I. The 
second constraint limits our approach to slow expansions of nuclear matter. Following Refs. 

see also the discussion in Ref. [jj on this point ), we assume that the relations 
fl23|) and for C^{t) and |-B k (t)| are valid also in unstable situations, provided that the 
conditions (3 < l/|r fc | and \Y k \ < vpk are still fulfilled, at least for time-averaged values. 
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The solution of the Fokker-Planck equation ()20|) is given by a Gaussian distribution. 
For simplicity we assume that the state of the system at t = is homogeneous on the 
average ( < Sg^t = 0) >= for k ^ ), in case of necessity a nonvanishing mean value 
could easily be introduced, then Eq. (|TH|) implies that this property holds also during the 
time evolution. By making the change of variables Sg^t) = Zy{t) / Ak{t) , we obtain for the 
probability distribution of density fluctuations the explicit expression 

F[*fik(t)] = #i(t)exp(- i£*£(f)i*gk(t)) , (27) 

with the variance cr^(t) given by 

n 2( f \ _ J1U _ nV 2 JoFk(t')dt> 2?r fc 1 /"* 2 f',r k (t")dt" T( . f) Mf) ,/ 

" fc(t) " ^ " 0)6 + ™ 2 F(fli) ^(*) Vo 6 T(t } (/"(f) - A(t')) ' 

(28) 

while iVi(i) in Eq. (J27)l is a time-dependent normalization factor. 

Equations (}2*T|) and (J28|l describe density fluctuations both in situations of local equilib- 
rium and for systems out of equilibrium. In the first case, for t ^> l/|Tfc| the variance relaxes 
towards the typical value of equilibrium thermodynamic fluctuations, with values of density 
and temperature appropriate to that instant, ( see Eq. (J24|) ), while in the second case 
the variance grows exponentially for fluctuations with wave number smaller than a certain 
time-dependent value. 

Starting from the probability distribution for density fluctuations given by Eq. (J27j) . we 
can determine the corresponding distribution for the size of the correlation domains. We 
would like to recall again that the stable and unstable cases can be treated within the same 
scheme, however here we investigate only the case in which the system has been quenched 
inside the spinodal zone. Then, it develops density fluctuations that grow with time and 
will eventually lead to the decomposition. 

The procedure followed in Ref. |l( for nuclear matter at constant density and temperature 
can directly be extended to the present expanding system, this is justified essentially by the 
Gaussian form of the probability distribution for density fluctuations. Then the distribution 
of the domain size is given by 



P(A,t)= — 



3v^F\ L(t) 



r °^) 4-5/6 p -[r (t)/L(t)]AV3 



where A is the number of nucleons contained in a correlation domain, L(t) is the length 
scale of the domain pattern and r (t) is the mean interparticle spacing. Like in Ref. 0], in 
order to take into account that A is a discrete variable we express the probability of finding 
a correlation domain containing A nucleons, Y(A,t), through the integral 

Y(A,t) = f dAP{A,t). (30) 

JA-l 

For large A, Y(A,t) tends to coincide with P(A,t). 



III. RESULTS 

A reasonable assumption for the initial condition is that the system can be modeled as 
nuclear matter at local equilibrium. This implies that the density fluctuations are those 
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FIG. 1: Time averaged growth-rate as a function of k at different times, from bottom to top 
t = 10, 30, 50 fm/c. The value of a is a = 0.01 

of the hydrodynamic regime. In particular the variance of fluctuations corresponds to an 
equilibrium state of given density and temperature at that instant, with erf (t = 0)\ eq given 
by Eq. (|24)1. Values of parameters T(0) = 6MeV and g(0) = 0.7g eq are compatible with 
the nuclear multifragmentation process. We also choose the initial expansion rate a = 
l'y(O) /t(0) I m the range ((1/6) 4 1CT 2 ) c/fm, which corresponds to a gentle expansion 
and cooling of the compound system formed after the collision. For a sphere of nuclear 
matter of radius R ~ 10 fm, the peripheral velocity lies in the interval (0.017 4- 0.1)c. Then 
we can investigate the transition from a quasistationary regime to a situation where the 
coupling between kinetic and hydrodynamic motions can give non-trivial effects. 

The time evolution of the fluctuations is ruled essentially by the growth rate r^.(t). This 
function depends on time through the product at, as a consequence the time average, 
Ff(t), scales with a. The magnitude of F% v (t) determines the range of validity of our 
approximations. The diffusion coefficients (J2rJ|) of the stochastic equation have been obtained 
in the hypothesis that thermodynamic fluctuations are the main effect, this requires r™(£) < 
T(t) and r™(t) < kvp- In Fig. we report the mean growth rate as a function of k at three 
different times for a given value of the parameter a ( a = 0.01 ). From that figure we can see 
that, for this value of a, the condition T^ v (t) < T(t) is satisfied up to a time t ^ 50 fm/c, 
moreover, the relevant values of the wave vector k and the values taken by the nuclear 
density in the interval < t < 50 fm/c are such that at each instant the quantity kvp is 
less than T. Thus the condition T™(t) < T(t) also implies T™(t) < kvp. According to the 
scaling law for the function T^ v (t), these limits are respected for longer times if a is smaller. 

In Fig. El the variance cr^(t) is displayed for three different values of the parameter a and 
at two different times for each value of a. The times have been chosen in such a way that 
the points representing the system on the phase diagram lie cose to each other inside the 
spinodal region in all three cases. The values of density and temperature at a given instant 
are lower for higher values of a. From Fig. El we can see that for slower expansions the 
variance is a faster function of time. Thus, at the same point inside the spinodal region, the 
density fluctuations are much larger for lower values of a. The slower the motion along an 
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FIG. 2: The variance as a function of k for three different values of a: solid lines, a = 10 -2 c/fm 
and t = 30,50fm/c; dashed lines, a = (1/3)10™ 2 c/fm and t = 60, lOOfm/c; dot-dashed line, 
a = (1/6)10~ 2 c/fm and t = 80, 150fm/c. The upper line corresponds to the later time. 



isoentrope, the easier is the decomposition of our nuclear system. 

The probability Y(A,t) of Eq. (|5Uj) is completely determined once the ratio between the 
length scale L(t) and the mean interparticle spacing ro(t) has been fixed. The parameter 
L(t) is the length characterizing the decrease of the correlation function with distance, while 
the correlation function itself is the Fourier transform of the variance cr\{t) with respect to 
k. A sufficiently accurate expression of the length L(t) is given by the simple relation 

m = M*) + A fc (t) ' (31) 

where /cm(^) is the value of k at the maximum of cr%(t) and A k (t) is the width of the broad 
peak in Fig. El This simple estimate of L(t) is in good agreement with the more precise 
evaluation made in Ref. |l| for the stationary case. By comparing Figs. ^ and El it can 
be seen that the positions of the maximum of T™(t) and of crf,(t) tend to coincide with 
increasing t. 

In Fig. El the ratio L(t)/ro(t) is shown as a function of time for three different values of 
the expansion rate a. The time interval for each of the three cases has been chosen so to 
allow for a comparison of the ratio L(t)/ro(t) for the same values of g(t) and T(t). In all 
cases the ratio L(t)/ro(t) approaches an almost constant value, lower than the initial one, for 
sufficiently long times. This fact can be ascribed to two competitive effects: with increasing 
t the broad peak of cr| (t) shown in Fig. Elgets narrower, at the same time &m(0 increases 
because the system reaches deeper regions inside the spinodal zone. The steeper initial slope 
seen in Fig. El instead, can be attributed to the fact that both /cm(^) and A k (t) increase with 
time, while cr^(t) changes from its initial value. This process is faster for rapid expansions. 
The ratio L(t)/r (t) is the only parameter contained in the probability distribution Y(A,t). 
The formation of the pattern of correlation domains can be viewed as a stationary process 
starting from the beginning of the plateau of L(t)/ro(t), even if the nuclear system is still 
expanding. In Fig. Elthe last point ( t = 70fm/c ) of the plot with a = 0.01 c/fm should be 
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FIG. 3: Behavior of the ratio between the length scale and the mean interparticle spacing L(t)/ro(t) 
for three different values of a: squares, a = 10~ 2 c/fm; circles, a = (1/3)10~ 2 c/fm; diamonds, 
a = (l/6)l(r 2 c/fm . 

taken with caution since it lies beyond the estimated range of validity of our calculations. 

In the present scenario, once expansion and cooling have brought the bulk of the system 
within the spinodal zone of the phase diagram, dynamical fluctuations grow with time until 
they cause the decomposition of the nuclear system. To be more definite, we fix the onset 
of the phase separation at the instant when the maximum of the variance cr|(t) becomes 
twice its initial value for k = 0. In addition, we identify the probability Y(A, t) with the size 
distribution of the domains containing the liquid phase. Our choice of the instant for the 
onset of the phase separation is also supported by the observation that this time roughly 
coincides with the beginning of the plateau observed in the behavior of the ratio L(t)/r (t). 
Thus, while fluctuations continue to grow with time, the size distribution Y(A,t) does not 
change considerably. The times for the beginning of the plateau in the three cases shown in 
Fig- El corresponding to a = (1, 1/3, 1/6)10~ 2 c/fm, are t ~ 50,80, 100 fm/c respectively. 

In order to assess the validity of our approach, we compare the results of our calculations 
with the corresponding experimental data by identifying the probability Y(A,t) with the 
distribution of the fragment yield observed in processes of nuclear multifragmentation. The 
probability Y(A, t) is calculated with the value of L(t)/r (t) at the plateau. Like in Ref. [l[, 
finite-size and surface effects are not included in our model, consequently a transition from 
surface to bulk multifragmentation, as suggested recently by the analysis of the ISiS results 
[l7| . cannot be described within our approach. However, we expect that our treatment 
should reliably account for the gross features of the bulk instabilities of nuclear matter, at 
least at a qualitative level. 

Since experimentally the fragments are detected according to their charge, we have to 
transform Y(A, t) into the corresponding function of Z, Y(Z, t). For this purpose, we assume 
a homogeneous distribution for the charge, Z = [(1 — I)/2]A, with I = (N — Z)/A and use 
/ = 0.2, that corresponds to the average asymmetry of the nuclear systems considered. 

In Ref. P it has been shown that the probability Y(Z) can be fitted with good accuracy 
by a power law Y(Z) = Y Z~ Te -t'f . The parameter r e ff depends only on the expansion rate 
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FIG. 4: Effective exponent of the power law, Y{Z) = YoZ~ Te ff, as a function of the collective 
radial energy per particle at the break-up (the corresponding range of the expansion parameter is 
(1/6)1(T 2 < a < 1(T 2 ). 

of nuclear matter, once the temperature and density at t — have been fixed. We have 
taken the initial values T(0) = 6 MeV and Qo(0) = 0.7 g eq , however, even if the temperature 
is changed in a rather large interval, the value of the ratio L(t)/ro(t) at the plateau does 
not change appreciably. In Fig. 0]the effective exponent r e ff is shown as a function of the 
radial collective energy per particle, Ehy, at the break-up of the nuclear system. This energy 
has a one-to-one correspondence with the parameter a, for given initial conditions, and is 
a more suitable observable than a for a comparison with the experimental data. Figure 0] 
shows a moderate rise of r e ff when the energy of the hydrodynamic motion increases. The 
cause of the expansion can be either dynamical, following an initial compression, or thermal, 
due to a release of energy in the bulk of the nuclear system from very energetic particles, 
for instance protons, antiprotons or pions. In both the cases the rapidity of the expansion 
and, as a consequence, the exponent r e ff increase with the excitation energy of the nuclear 
system. 

In Fig. 03a comparison is made between the charge distributions predicted by our approach 
and recent experimental data obtained by the the INDRA Collaboration for l29 Xe + Sn 
and 155 Gd + U collisions at E = 32 A MeV and at E = 36 A MeV respectively [pi . In 
Fig. the same comparison is made for Xe + Sn collisions at E = 50 A MeV |l4j. Our 
calculations have been performed for three different values of the parameter L(t)/r (t), the 
corresponding values of the effective exponent r e ff are given in the legend of Fig. |SJ We 
have normalized both the experimental and the calculated distributions to one in order to 
perform the comparison on an absolute scale. The agreement between experimental data 
and our calculated charge distributions is quite satisfactory for Z < 30 4- 35 at the lower 
incident energies of Fig. and for Z < 20 -f- 24 at the higher energy of Fig. |BJ For the higher 
values of Z the observed distributions display a steeper slope than that predicted by our 
calculations. This faster decrease of the experimental yields should be ascribed to finite-size 
effects |22| that have not been included in our nuclear matter treatment. The comparison 
also show that experimental data at higher energy are better reproduced with smaller values 
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FIG. 5: Comparison of fragment distribution Y(Z) calculated for L(t)/ro(t) = 1.3 ( dot-dashed 
line ) 1.6, ( dashed line ) and 2 ( solid line ) ( r e ff = 1.65, 1.51, 1.41 respectively ), with experimental 
distributions for the reactions 129 Xe + Sn at E = 32 AMeV (triangles) and 155 Gd + U at E = 
36 AMeV (circles) ^f. Both data and curves have been normalized to 1. 
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FIG. 6: Same of Fig. El but for the reaction Xe + Sn at E = 50 AMeV (diamonds) Q. 

of the parameter L(t)/ro(t) and correspondingly larger values of T e ff, compared with the 
values needed to reproduce the data at lower energy. This trend is explained by our scheme, 
since a growth of the radial collective energy corresponds to an increase of r e // ( see Fig. El ) 
and it is reasonable to expect that for central collisions the radial collective motion becomes 
faster when the incident energy increases. 

A further remark about the effective exponent T e ff. It has been experimentally observed 
that the onset of an expansion in nuclear matter is a signal for a transition from surface- 
dominated to bulk emission, as expected for a spinodal decomposition, and that the power- 
law parameter r e jf grows with excitation energy for expanding nuclear matter [2, [J. Our 
results are in qualitative agreement with this observed behavior and so is our estimate of 
the radial collective energy. However, the experimental values of r e ff reported in Ref. 
are systematically larger than our values. 
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IV. SUMMARY AND CONCLUSIONS 



We have studied the density fluctuations of expanding nuclear matter within a one-body 
treatment of nuclear dynamics. In our approach the fluctuations are generated by adding a 
stochastic term to the mean field and they develop about a local-equlibrium state described 
by a hydrodynamic expansion. This additional random force is constrained by an extension 
of the usual fluctuation-dissipation theorem to the case of uniformly expanding systems. In 
this more general context, we have confirmed the main results obtained in Ref. about 
the particular physical conditions in which a white-noise assumption for the stochastic field 
can be retained: the density and temperature of the system should be sufficiently close to 
the borders of the spinodal region in the (q,T) plane, so that the limit u/T <C 1 gives 
a reasonable approximation to the density-density response. In this case the equilibrium 
fluctuations can be adequately described by means of thermodynamic functions. 

We have extended the results obtained for the probability distribution of fluctuations in 
a system at local-equilibrium to unstable situations. This has been done by extrapolating 
the relevant quantities across the boundary of the spinodal region. Because of the linear 
approximation used for evaluating the response of the system to the stochastic force, the 
probability distribution of the fluctuations is Gaussian. The variance of the distribution at a 
given point in the (g, T) plane, strongly depends on the expansion velocity and instabilities 
are favoured in the case of slower expansions. 

In Sec. Ill Bl we have discussed a procedure to determine the size distribution of the 
domains containing correlated density fluctuations. This procedure is quite general and 
can be applied to any Gaussian distribution. For the size distribution we have found an 



expression similar to that obtained in Ref. It again contains only one parameter, 

the ratio between the dynamical correlation length and the mean interparticle spacing, 
L(t)/ro(t). However, in the present case the time behavior of L(t) is mainly determined by 
the expansion velocity, in other words by the collective energy stored in the excited nuclear 
system. In particular, the parameter L(t)/ro(t) displays a flat behavior in time after a 
certain instant that depends on the expansion velocity. Roughly at this time, the amplitude 
of fluctuations acquires at least twice its initial value and keeps growing. Therefore, the 
decomposition of the system can be viewed as a quasi stationary process. 

In Sec. Illll we have compared the obtained mass distribution to the yield of intermediate- 
mass fragments observed in the multifragmentation of heavy nuclei. The comparison with 
experiment has shown that our approach fairly reproduces the measured charge distributions 
with Z < 30 -T- 35 for collisions at lower energy and with Z < 20 at higher energies. Since we 
have modeled the system as infinite nuclear matter, we expect to overestimate the number 
of large fragments. However, it could be worthwhile studying whether and how the finite 
size effects could become more important with increasing incident energy. 

Our approach does account for both the power-law distribution found experimentally 
and the observed behavior of the effective exponent with the excitation energy, but the 
calculated exponent is systematically smaller than that determined from experiments. This 
discrepancy deserves further investigation. 

APPENDIX 

In this Appendix we briefly discuss the statistical properties of nuclear matter that un- 
dergoes the expansion described in Sec. Ill Al For this hydrodynamic motion we can assume 
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that in the local reference frame moving with the hydrodynamic velocity u(r, t), the physical 
system is described by the equilibrium density matrix: 

m = W) e ~ mH ' (A,1) 

where Z(t) = Trg(t) is the partition function and H is the Hamiltonian in the local reference 

fl*ctIXXG 

f?=£ (p,-mu(r J -, i y + l £v , (rj ^ ) (A2) 

It can be shown that, for a generic couple of states \m > and |n > 

< m\H\n >=< m\ e^*^ H e^^^n > , 

where x( r , t) is a gauge function such that Vx(r, t) = mu(r, t). Thus, the statistical density 
matrix of Eq. (|A.1)1 is equivalent to 

£(t) = ^ e ffl i x(r i ,t) e -^P e -iE jX (r j ,t) _ (^3) 

z(t) 

Because of the properties of the trace the partition function 

Z{t) =Tre-^ H 

is independent of the velocity field. The density matrix of Eq. (|A.3|) reproduces the correct 
expression both for the collective current and for the density of the hydrodynamic kinetic 
energy 

\ < E (- S(t - r,) + 6{r - r,) El) >= g ( t )u(r, t) 

and 

< E 5 ( r - r i)) >= 7^oWu 2 (r, t) + , 

where Skin is the density of the intrinsic kinetic energy. 

The observables that commutate with r take the same mean value as for a homogeneous 
medium at equilibrium, having the specific values of density and temperature at a given 
instant. The same result holds for the equal-time correlation functions of such observables. 
In particular, for the density fluctuations we obtain 

< 6q(t, t)5g(v', t)>= -^-Tr Se^Sgifje-^) 11 . (A.4) 

As far as the correlation functions at different times are concerned, the result of Eq. (jllj) 
of Sec. Ill Al suggests that, when the collective expansion is slow with respect to the Fermi 
velocity, these functions depend only on the values of density and temperature of the system 
at a given instant. In this case, we can reasonably neglect the velocity field J2j x( r j; t) in the 
expression of the density matrix when evaluating the density correlations. Then we expect 
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that a instantaneous fluctuation-dissipation relation still holds for a medium undergoing a 
hydrodynamic uniform expansion. In fact, it can be shown that 

2 

< (tfgktf0-k)H >f= ~ 1 _ e _ /?(f> Im£> fc (o;,f) , (A.5) 

where the time Fourier transforms have been evaluated respect to t — t'. The statistical 
average is calculated with the density matrix 

g[t) Z(t>) 
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